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Abstract 

Electronic structures of the zigzag bilayer graphite nanoribbons(Z-BGNR) with various ribbon 
width N are studied within the tight binding approximation. Neglecting the inter-layer hopping 
amplitude 74, which is an order of magnitude smaller than the other inter-layer hopping parameters 
71 and 73, there exist two fixed Fermi points ±k* independent of the ribbon width with the peculiar 
energy dispersion near k* as e(k) ~ ±(k — k*) N . By investigating the edge states of the Z-BGNR, 
we notice that the trigonal warping of the bilayer graphene sheets are reflected on in the edge state 
structure. With the inclusion of 74, the above two Fermi points are not fixed, but drift toward 
the vicinity of the Dirac point with the increase of the width TV as shown by the finite scaling 
method and the peculiar dispersions change to the parabolic ones. The edge magnetism of the 
Z-BGNR is also examined by solving the half-filled Hubbard Hamiltonian for the ribbon using the 
Hartree-Fock approximation. We have shown that within the same side of the edges, the edge spins 
are aligned ferromagnetically for the experimentally relevant set of parameters. 

PACS numbers: 73. 21. Ac, 73.9Q.tf, 73.21. -b 
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I. INTRODUCTION 



Graphite nanoribbons(GNR) are the quasi one-dimensional carbon sheets, whose width 
and the chirality can be well controlled by the current nano-lithography technique. [l| Due to 
the high mobility of the sample and the relativistic energy dispersion relation of the massless 



Dirac fermion, there has been a lot of interest both theoretically and experimental 



tigate the electronic properties of the GNR. 
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The electronic properties of the GNR are closely controlled by the chirality as is well known 
for the case of carbon nanotube. Among the GNRs with different chiralities, the armchair 



and zig zag GNR have been studied most intensively for the decades. B,y,Q,y,y,y, 
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171 | Unlikely from the carbon nanotube with the compact angle variable 



for the transverse modes, the GNR can support the edge states. Recently, the edge state of 
zigzag GNR has been of great interest due to its peculiar dispersion relation with the almost 
flat edge bands near the zero energy. 



.ay , [ijj] When the Coulomb interactions are taken 



into account, the existence o 
kinds of ribbon edges 



121 Q, 



;he flat edge bands may lead to the edge magnetism for various 
| They have shown that the edge magnetism of zigzag 
GNR is ferromagnetic along each edge while anti-ferromagnetic between the two different 
edges due to the bipartite nature of the lattice structure. This kind of edge magnetism has 



been widely studied within first-principle calculations. [20j, 



21 



22 



23] Furthermore, its po- 



tential as a half-metallic material under lateral electric field has been demonstrated within 



a first principle approach by Son et. al. [18] Hence to identify the edge states of the various 



GNRs is crucial to understand both the magnetic and electric properties of the system. 

In the paper, we study the electronic and magnetic structure of the zigzag bilayer graphite 
nanoribbons with various ribbon width N. First, we give a brief review of the bilayer 
graphene system. The bilayer graphene system is two coupled graphene sheets and its 
tight binding electronic structures are determined by four leading characteristic hopping 
parameters as shown in Fig. 1(a). It has demonstrated quite unique electronic properties 
compared to the single graphene such as the larger Hall step and different kinds of chiral 
quasiparticles and the distinct electronic structures. [24J The most distinct electronic feature 
of bilayer graphene system is the very existence of the trigonal warping of the Fermi level 
in the vicinity of the Dirac points as demonstrated in the inset of Fig. 2(a). It means that 
the band of the graphite bilayer near the Dirac point does not form a cone but a trigonally 
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curved cone and there exist three additional pockets along the directions of zigzag chirality. 
It has been shown that these peculiar electronic structures, which are controlled by the 
inter-layer hopping parameters 71 and 73, produce extra degeneracies in the Landau levels. 
The three extra dirac cones are known to play an important role in determining the minimal 
conductivity of the graphite bilayer. 25j 

In the work, we focus on the edge state structure of the Z-BGNRs with various ribbon 
width N. Since the three legs in the band structure of the graphite bilayer lie along the di- 
rection of zigzag chirality, we have chosen the Z-BGNR in order to see the possible signature 
of the trigonal warping imprinted on the edge state structure of the system. The schematic 
structure of the Z-BGNR is shown in Fig. 1(a). We have chosen the Bernal stacked graphite 
nanoribbons, where the widths of the two layers are the same and the width is defined to 
be the number of dimer lines N as shown in Fig. 1(b). We use the tight binding method 
to calculate the band structures of the Z-BGNRs with varying widths. According to the 
experimental data, the inter-layer hopping amplitude 74 is much smaller than 73 by an 
order of magnitude. j26j Neglecting 74, we notice that the Z-BGNR has two fixed Fermi 
points at k = ±k* with k* = 2 cos" 1 (\/7i7}/2) and the following energy dispersion relation 
e ~ ±(k — k*) . As N goes to infinity, the dispersionless edge band appears. By analyzing 
the semi-infinite Z-BGNR case analytically, we have demonstrated that the trigonal warp- 
ing manifests itself in the edge band structure of the Z-BGNR. Interestingly, the existence 
of the edge states seems to be consistent with the condition imposed by the Berry phase 
argument. With the inclusion of 74, the particle-hole symmetry is broken and the energy 
gap appears in the dispersion. In this case, the conduction band minimum and the valence 
band maximum occurs at two different points k e and k h respectively. They are no more fixed 
points independent of iV rather drift toward Dirac point with increasing N. By examining 
a semi-infinite Z-BGNR with 74 taken perturbatively and using the scaling method as well, 
we have analyzed the edge states in further detail. We notice that one edge band remains to 
be flat while the other edge band becomes dispersive. These two edge bands appear below 
e = and split proportional to the hopping parameter 74. Finally, the magnetic properties 
of the Z-BGNR are studied by investigating the Hubbard Hamiltonian obtained from the 
inclusion of the on-site Coulomb repulsion U based on the Hartree-Fock approximation. We 
have calculated the magnetic moments at the edge as a function of U for half-filled system 
at zero temperature. Our main concern is that the possible half metallicity would persist in 
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the Z-BGNR upon lateral gate bias. Our result has shown that for the realistic value of U 
(C//70 — 1), the ferromagnetic alignment of the edge spins within the same side of the edges 
is favorable. This implies that the spintronic application of the Z-BGNR is quite promising. 

The paper is organized as follows. In Sec. II, a detailed descriptions about the tight binding 
model for the Z-BGNRs are given. In Sec. Ill, we investigate the electronic structures of Z- 
BGNRs with or without the small inter-layer hopping parameter 74. We focus on the edge 
state structure based on the formulae of Sec. II. Without 74, the coupled equation to obtain 
the energy dispersion relation can be simplified and one can extract the important features of 
the edge states in the system. We analyze the semi-infinite Z-BGNRs to obtain the analytic 
formulae of the edge states amplitudes. In Sec. IV, the interlayer edge magnetism of the 
Z-BGNR is studied by including the Hubbard type onsite-interactions. Finally, summaries 
will be given in Sec.V. 



II. TIGHT BINDING HAMILTONIAN 



We introduce the tight binding model for the Z-BGNR by considering the four hopping 
parameters 70,71,73 and 74, where 70 is the intra-layer hopping and the others the inter- 
layer hopping amplitudes as shown in Fig. 1(a). This type of approach has been successfully 
employed to study the electronic structure due to the p^-orbital of various carbon based 
materials including the GNR. {3], Q] We begin with the following Hamiltonian 

H = -J2u J c\c j (1) 

with 

t tj = {<p{r f -R<)\H\<p{?-R j )) (2) 

where the summations over i and j include all the hopping processes related to the four 
hopping parameters. For the sake of simplicity, we redefine the index i in terms of the 
following three indices such as q = c a (m, A), where a represents the unit cell of Z-BGNR, 
m the index of dimer lines which runs from 1 to N, and A one of the four sublattices A,B,A 
and B as shown in Fig. 1(b). By transforming this to the collective modes 

c m)X (k) = 4^Y1 e"^ Q c Q (m, A) (3) 

* a 



where y a denotes the position of a-th unit cell. The general states can be written by 



By substituting Eq. (J3J) into Eq. (pQ), one can obtain the following four recurrence relations 



eip m A = 7o(2 cos -Vwb + ipm-w) + 7s(2 cos -V> m _iB + ^Vb) + 74(2 cos -ip mA + ^ m _ lA f>) 



e ^mA = 7o(2 COS -1p mS + Vm-lfi) + ll^mB + 74(2 COS -lfj mA + </Wia) 

k /c Aj 

e ^mB = 7o(2 COS -iJ mA + ^ m+lA ) + 73(2 COS -^m+lA + Vw) + 74(2 COS -^mB + VWlZ?)- 

Here we use the dimensionless wave vector k in which the primitive lattice vector of 
the ribbon is embedded so that the range of k runs from —it to ir. Since the Z-BGNR is 
terminated along the x-axis, the following boundary conditions are imposed: ip \ = ipN+i\ = 
0. From the recurrence relations, we notice that, when 74 is taken to be zero, particle-hole 
symmetry exists in the band structures. In other words, for 74 = 0, one set of eigenvalue 
and eigenvector solution of {e, ip m A, ^mB-, ^Vb} guarantees the existence of the other 
orthogonal solution {— s, —ipmAyipmB, — ^mA'^/Vsl- Hence the band structure is symmetric 
with respect to e — 0. On the other hand, when 74 is included, the particle-hole symmetry 
is broken. 

III. ELECTRONIC STRUCTURE OF THE Z-BGNR 

We study the electronic structures of the Z-BGNRs with finite ribbon width iV with spe- 
cial focus on the edge states. We first investigate the semi-infinite Z-BGNR and analytically 
show the signatures of the trigonal warping of the Fermi surface in the edge states. Hence 
our main interest lies in the bands near e = 0, which contain the states localized at the 
ribbon edge. For e = 0, the four coupled equations in Eq. (5) are decomposed into the two 
coupled ones. One of the two coupled equations only involves ipB and ipg. By solving these 
equations, we notice that they lead to unphysical solutions unless ip m B = Vw? = 0- The 




(4) 



£i>mB 




k 
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other two equations, which contain tpA and ip A , are written by 

i>m+2A +P(2- 7i7 3 )^m+lA + (P* - lll^mA = (6) 

where 71 = 71/70, 73 = 73/70) an d /3 = 2cos(A;/2). The first equation can be solved by 
substituting an ansatz of ip m A oc /i m into Eq. (6). As expected from the double degeneracy 
of band structures of bilayer graphite at the particle-hole symmetry (e — 0), one can obtain 
two independent solutions /i± as follows 

„ ± = -0(1 - 2^£) (7) 

±y / fif 3 (i-/? 2 ) + /J 2 (^) 2 . 

where the localization lengths are given by £± = | In |//±|| _1 . 

We denote the two solutions as for /x± respectively. In order to have a physical 
normalizable solution, one needs to impose the following condition for /i: < 1. This 
condition restricts the range of the k values supporting the zero energy edge states as follows 

For ^ mA , 2 cos- 1 V1+ o flf3 < k < vr (8) 



For ^ 



mA 




2 < k < 7T. 



At fc = 7r, the edge state has a form ip mA = (± v /7 1 7 3 ) m and = — (±\/7T7^) m+1 /7i- 
Hence the localization length of the semi-infinite Z-BGNR is finite at the zone boundary, 
while the semi-infinite monolayer GNR is totally localized at one edge for k = 71. It is 
interesting to notice that the Z-BGNRs also support a completely localized mode at k* = 
2 cos" 1 v /7i7 3 /2 at which /i — 0. 

In Fig. 2(b), the schematic diagram of the edge states are plotted with respect to k, 
where the two edge modes are represented by two different colors. The yellow(pale gray) 
edge mode V'mA ex ists from the point 7\ at k — 2 cos -1 y/l + 7173/2 to the zone boundary as 
observed in the monolayer GNR as well. Remarkably, we notice that the red(dark gray) edge 
mode ip^A h &s a forbidden region between the Dirac point D at k = 2n/3 and the point T 2 
at k = 2 cos _1 (l — 7i7 3 )/2. The existence of the forbidden region is the unique feature of the 
Z-BGNR. We draw two symmetric warped bands within this region, since the bands should 
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be continuous and particle-hole symmetric. By comparing to the graphene bilayer system, 
the point D corresponds to the Dirac point while the points T\ and T2 correspond to the 
three legs' Fermi points Li,L 2 and L 3 of the graphene bilayer which are indicated in the inset 
of Fig. 2(a).[24j| This correspondence can be explicitly shown by comparing these points. 
The previous study on the graphene bilayer has shown that the three legs are displaced by 
27173/ v^3 from the Dirac point D and are placed at the points with triangular symmetry. If 
one projects these leg points onto the axis of zigzag chirality, the two leg points are located 
at k = 2tc/3 + 27i7s/v / 3 and k = 2n/3 — 7173/ v3. These values previously obtained by the 
k ■ p approximation perfectly agree with the two points T\ and T2 in Eq. (8) for small 71 
and 73. We have confirmed the above feature of the zero energy edge states and the warped 
bands by numerically solving Eq. (5) for relatively wide Z-BGNR with widths N = 100 
and 200 as shown in the inset of Fig. 2. In order to see the above feature more clearly, we 
have taken about five times larger values for 7i and 73 than the experimental ones. Our 
numerical results agree well with that of the semi-infinite Z-BGNR. 

Recently, Ryu et. al. suggested a criterion for the existence of the edge states based on 
the bulk energy dispersion, which may be applicable to a certain class of edge states 27]. 
They have studied the following particle-hole symmetric Hamiltonian 

H = E = E 4(^) • ( ) 

k k 

where ct = (ct , c^) and a is a three component vector of the Pauli matrices. This type of 
Hamiltonian includes the Bogoliubov-de Gennes and the GNR Hamiltonian as well. Their 
criterion for this Hamiltonian to have zero-energy edge sates is that the closed trajectory of 
R(k) being confined in a 2D plane encloses the gap-closing point. We want to understand 
our result of Z-BGNR based on the above criterion. By analyzing the 2D graphite bilayer 
system using exact diagonalization method, one can obtain four energy bands, which we take 



to be |it!(fc)|.[27j] We have analyzed the behavior of one of |i?(/c)|'s and noticed that it goes 
to zero three times at k y = 2cos~ 1 y/\ + 7173/2, 2tc/3, and 2cos~ 1 (l — 7i73)/2 respectively. 
We presume that these three points indicate the boundaries of the three distinct topological 
sectors of |-R(&;)| where the closed loops within 27r/3 < k y < 2cos~ 1 (l — 7i73)/2 exclude the 
gap-closing point. Hence we argue that it corresponds to the edge mode of t/>~^. Strictly 
speaking, our bilayer graphene system is out of their scope, since the bilayer Hamiltonian is 
described by the 4x4 traceless hermitian matrix. Furthermore, since only the magnitude 



of R(k) are available, we are not able to directly confirm our conjecture yet. 

We now calculate the full electronic structures of the Z-BGNRs with finite ribbon width 
iV by numerically solving Eq. (5). We obtain the band structures of two finite width Z- 
BGNRs with N = 4, 5 as shown in Fig. 3. In order to study the edge states and the low-lying 
energy excitations analytically, we diagonalize the 4N x 4N matrix of Eq. (5) with e — 0. 
By taking the determinant of the matrix to be zero, that is, [(2 cos(A;/2)) 2 — 7i73] 2Ar = 0, 
one can see that the band touches e = at ±k* independent of N. In order to calculate 
the dispersion relations around (k,e) = (k*,0), we expand the determinant with respect to 
(k*, 0) and obtain the following results 



± v /^(2(* -*•))" (iV = odd) 



±y \{2{k - k*)) N (iV = even) 

where e is written in units of 70 . These peculiar dispersions are demonstrated in our nu- 
merical result as shown in the insets of Fig. 3. When the ribbon width iV goes to infinity, 
the energy dispersion becomes flattened leading to the localized edge states. The similar 
behavior has been observed for the monolayer GNR with k* being replaced with k — ir. 
Unlikely from the monolayer GNR, a small gap opens at the zone boundary. The magnitude 
of the gap can be calculated approximately by considering the weak coupling between those 
edge modes localized at the different edges. We define the following four states which are 
constructed from the two semi-infinite Z-BGNRs localized at the opposite edges for k = n 

\L) = ^ 1A (|iA)-7 3 |2i)+7 1 7 3 |3A) + ---) (11) 
\L) = ^(|li)-7 1 |2A)+7 1 7 3 |3i) + ...) 
\R) = if> NB (\NB)-7 3 \N-lB) 

+ia 3 \N-2B) + •••) 
\R) = ^(|L4>-7 3 |iV-LB> 

+7i73|iV- 2B) + •••), 

where \L) and \L) represent the localized states at the left edge, while \R) and \R) are those 
at the right edge. Using these four orthogonal states as the basis, the matrix elements of 
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the tight binding Hamiltonian can be calculated leading to the following results 

±(7173)^" (N = even) 



__ __ ( 12 ) 

±71(7173) 2 , ±73(7173) 2 (N = odd). 



N-l 



For the even width ribbons, the two split edge bands in the first quadrant converge to a 
single value e = (TlTs)^ 2 a t the zone boundary, while they remain to be split for the odd 
ones. This is due to the reduced symmetry of the odd width ribbons compared to the even 
ones. 



According to the previous experiments, the band parameter 7, 



24 



23, 



23, 



30 



was assumed to be much 



31| The typical values for 



smaller than 71 and 73 by an order of magnitude, 
the hopping amplitudes are taken to be 7i — 0.12, 73 — 0.1 and 74 — 0.014. [26| Hence it 
is of practice to neglect the contribution of 74. Here we want to investigate the effect of 74 
on the low-lying energy dispersion. Within the k ■ p approximation, we notice that not only 
the Fermi energy e = has a downward shift by Ae ~ — 27i7374> but also the particle-hole 
symmetry is broken upon including 74. 

Figure 4. (a) and (b) show the band structures for iV = 4, 5 obtained by numerically 
solving Eq. (5). One can observe the downward shift of the Fermi energy as expected. 
Interestingly, the degenerate two bands at k = k* for 74 = are split into the two bands 
with an indirect energy gap. We denote the extremal points as k e for electron pocket and 
k h for hole pocket as indicated in the insets of Fig. 4. When N goes to infinity, the bands 
become flattened. 

We have closely examined the positions of k e and k h with increasing N and observed that 
the inclusion of 74 leads to a qualitatively different behavior. In Fig. 5, we plot k e and k h 
with respect to 1/N and extrapolate to iV — > 00. With increasing N, k e and k h , which used 
to be fixed at k*, are moving towards the vicinity of the Dirac point at k = 2n /3. By fitting 
the energy dispersions near k* e and k* h up to iV = 4, we find that the peculiar dispersions of 
e ~ (k — k*) N reduce to the parabolic one except for the N = 1 case. 

IV. MAGNETIC PROPERTIES OF THE Z-BGNR 

When the Coulomb interaction U is present in the monolayer graphene system, sublattice 
ferromagnetism can appear beyond a certain critical value of U. The most favorable spin 
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configuration is known to be the opposite magnetization between two sublattices A and B 



which is consistent with Lieb's theorem regarding the magnetism of bipartite lattices. |32| 
For the case of the Z-GNR, the existence of almost flat bands at zero-energy stabilizes the 
sublattice edge ferromagnetism even for the weak Coulomb interaction. Recently it has 
been proposed that the Z-GNR can become half-metallic when the transverse electric field 



is applied. [17|, [18|, |22j We have studied the edge magnetism for Z-BGNR by applying the 
Hartree-Fock approximation to the half-filled Hubbard Hamiltonian. Following the Hartree- 
Fock decoupling, the Hubbard Hamiltonian can be written by 



H M f = ~ l A c i° + U Yl 



(13) 



where we have included the on-site Coulomb interaction U and n ia = c\ a Ci tCT . The mean 
electron density (rij CT ) is kept to be constant in a given dimer line due to the translational 
invariance. Based on the above Hartree-Fock Hamiltonian, we have solved the self-consistent 
equations. Since the inter-layer hopping parameters 71 and 73 are much smaller than the 
intra-layer one 70 by an order of magnitude, we take the initial configuration for itera- 
tion as the two coupled ground state configurations of the monolayer GNR. Following the 
symmetry of Z-BGNR, there exist two possible inter-layer spin configurations as shown in 
Fig. 6. One is the anti- ferromagnetic alignment (AFM) of inter- layer edge spins within the 
same side of edges: {n A , m ,a) = (™s,7v- m +i, CT ) and (n B , m ,a) = (n^ N _ m+x ^). The other is 
the ferromagnetic alignment (FM) of inter-layer edge spins: (nA,m,a) = (n§ jv-m+i s) anc ^ 

(n>B,m,o) = \ n A,N-m+l,a) ■ 

In order to see which configuration is energetically favorable, we have calculated the total 
energy of the system in those two cases. The differences between two energies Eafm — Efm 
are plotted as a function of U in Fig. 7. Here we have examined a Z-BGNR of width 
N = 10 and the inter- layer hopping parameters 71 = 73 = 0.1. We notice that both AFM 
and FM configurations open a gap by lifting off the flat bands. Figure 7 demonstrates that 
the FM configuration becomes immediately stable upon including the Coulomb interaction. 
Since the ferromagnetic alignment of the upper and lower edge spins leads to the inter- 
layer anti-ferromagnetic coupling between the nearest neighbor spin pairs connected by 71 
and 73, the ferromagnetic alignment will become energetically favorable for large U. Hence 
for the ex per imentally relevant value of U, that is, U/'yo ~ 1, the FM configurations are 



stable. 



33j, 



34 



35| We have also studied the bulk magnetism for the bilayer graphene system. 
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In the inset of figure 7, the energy differences Eafm — E FM are plotted as a function of U. 
With the increase of U above U c = 2. 570, the paramagnetic spin state becomes the FM one, 
which agrees with the result of the monolayer graphene system. Hence for the Z-BGNR, the 
spin fluctuations at the edges enhance the appearance of the FM configurations even below 
U c . 

V. CONCLUSIONS 

We have studied the electronic and magnetic structures of edge states of the zigzag bilayer 
graphite nanoribbons within the tight binding approximation. Neglecting the band parame- 
ter 74, we obtained a peculiar dispersion relation e ~ ±(k — k*) N with k* = 2 cos" 1 ^/^%/2 
which is fixed independent of the ribbon width N. As iV increases, the dispersionless edge 
bands at e = will appear just like a monolayer zigzag GNR. By investigating the semi- 
infinite Z-BGNR analytically, we notice that the trigonal warping of the bilayer graphene 
leads to the following interesting effect on the edge state structure of the Z-BGNR. One 
of the edge state at e = is absent within the region between Dirac point (A; = 2ti/3) to 
k — 2 cos" 1 (l — 7173/2), where the bulk warped bands used to be present. With the inclusion 
of 74, the particle-hole symmetry is broken in our Z-BGNR system. There exists an indirect 
gap between two edge bands, whose band extrema are located at k e and k h respectively. 
Instead of the fixed k* as for the case of 74 = 0, k e and k h approaches near to the Dirac 
point with increasing N. The magnetic property of the Z-BGNR is studied by solving the 
Hubbard type Hamiltonian using Hartree-Fock approximation. We have shown that for the 
realistic value of U, the ferromagnetic alignment of inter-layer edge spins within the same 
side of edges is favored. 
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FIG. 1: (Color online) (a) Our model system of Bernal stacked Z-BGNR. Both the three inter- 
layer(7i, 73 and 74) and one intra-layer hopping(7o) parameters are shown. The Z-BGNR has a 
translational symmetry along T. The sublattice A is colored yellow(pale gray) while B is colored 
green(dark gray), (b) The definition of the width of the Z-BGNR. Here the red lines(pale gray) 
are the lower graphene sheet and the black lines the upper one. The tilde above A and B stand 
for the upper layer and the black box represents a unit cell of the Z-BGNR. 

FIG. 2: (Color online) (a) Projected band structure of 2D graphite bilayer along the direction of 
zigzag axis. The boxed region near the Dirac point is magnified in the inset, which exhibits the 
trigonal warping. D is the Dirac point at k = 2ir/3 and three L points are the Fermi points of 
three nearby pockets, (b) A schematic diagram of the zero-energy edge states of the semi-infinite 
Z-BGNR near the Dirac point within the red box in Fig. 2(a). At e = 0, the two eigenstates are 
drawn with different colors. Red(dark gray) for and yellow(pale gray) for D is the Dirac 
point while T\, T2, and the two warped bands reflect on the effect of the trigonal warping of the 
graphite bilayer. The inset shows the energy dispersion curve obtained by numerical method for 
the Z-BGNRs with finite widths N = 100, 200. 

FIG. 3: (Color online) (a) The band diagram of the Z-BGNR with the following parameters 71 = 
0.12, 73 = 0.1 and 74 = for TV = 4. The inset, a magnified view of the low energy region, shows 
the power law dispersion of the edge bands, (b) The band diagram of the Z-BGNR for N = 5. 

FIG. 4: (Color online) (a)The band diagram for the Z-BGNR with the band parameters 71 = 0.12, 
73 = 0.1 and 74 = 0.014 for N = 4. The inset, a magnified view of the low energy region, 
demonstrates the broken particle-hole symmetry and the parabolic dispersions, (b) The band 
diagram of the Z-BGNR for N = 5. 
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FIG. 5: (Color online) The plot of k e and k with respect to l/N. The extrapolated values of k e 
and k h approach near to the Dirac point (k = 2ir/3 ~ 2.094). 

FIG. 6: (Color online) Two possible configurations of the magnetic moments are illustrated at the 
left part of the Z-BGNR. Black and gray lines represent dimer lines of the upper and lower layer 
each. Yellow(pale gray) arrows are the magnetic moments of the lower layer while red(dark gray) 
arrows for upper one. (a) The AFM configuration. At the left edge, the spin polarization at upper 
edge(lyl) and lower edge(lyl) layer are opposite, (b) The FM configuration. The spin polarization 
at upper edge(1^4) and lower edge(1^4) layer are parallel. 

FIG. 7: (Color online) A plot of the energy differences per unit cell between the AFM and FM 
configurations for the N=10 Z-BNGR and 2D graphite bilayer( inset) plotted as a function of U. 
The unit of the energy is 70. 
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